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ANALYSIS OP FULLY DEVELOPED TURBULENT HEAT TRANSFER AT LOW PECLET NUMBERS 
IN SMOOTH TUBES WITH APPLICATION TO LIQUID METALS 
By Robert G. Deissler 


SUMMARY 

An analysis was made of beat transfer for fluids flowing turbulently 
at low Peclet numbers in smooth tubes . Previous analyses for flow of 
gases and. liquid, metals at low Peclet numbers gave higher heat-transfer 
coefficients than were indicated, experimentally. When the mixing-length 
theory was modified, in order to account for the heat transferred by 
conduction to a turbulent particle as it moves transversely, the pre- 
dicted results were brought into agreement with the experimental results. 


INTRODUCTION 

Analyses of heat transfer for liquid metals flowing turbulently in 
smooth- tubes (lqw Prandtl numbers) are given in references 1 and 2. The 
predicted heat-transfer coefficients from these analyses are considerably 
higher than those determined in the experimental heat-transfer inves- 
tigations for mercury -and lead, bismuth given in references 2 to 4. It 
should be noted, however, that the experimental data were all obtained at 
low Peclet numbers. If the analytical and experimental results for flow 
of gases are compared at low Peclet numbers, they are also found to be 
in disagreement (references 5 and 6) . These results indicate that the 
eddy diffusivities for momentum and heat transfer cannot be considered 
equal at low values of Peclet .number . The fact that the eddy diffusivities 
for heat and momentum transfer cannot be considered equal at low Peclet 
numbers is also illustrated by the temperature distributions given in 
reference 5 where a considerable separation of data was noted at low 
Peclet numbers. Analytical and experimental investigations of the 
variation of the ratio of eddy diffusivities for heat and momentum transfer 
are given in references 7 and 8. In an analysis made at the NACA Lewis 
laboratory, the heat transferred by conduction to a turbulent particle 
as it moves transversely Is considered so that the eddy diffusivity for 
heat transfer is reduced to a value below that for momentum transfer. 
Although a simplified theory is used in the analysis, it appears to be 
sufficiently accurate for predicting heat transfer at low Peclet numbers . 
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ANALYSIS 

When turbulent flow exists in a tube, portions of the fluid move 
about in random fashion. If heat transfer between the tube wall and 
the fluid tales place, a temperature gradient occurs across the tube 
and some of the portions of the fluid move transversely into regions of 
different mean temperature . This motion produces heat transfer in 
addition to the transfer that tales place by local molecular conduction. 

In the usual mixing -length theory of turbulent heat transfer, an 
eddy or turbulent particle is assumed to originate from an instability to 

at a point in the fluid having the mean temperature of the fluid at the cG 

point, to travel , transversely at constant temperature a distance of one 
mixing length, and then to mix with the fluid and assume the average 
temperature of the fluid at the point of mixing (reference 9) . (In 
the present discussion, an eddy or turbulent particle is considered to 
be a portion of fluid all parts of which move with approximately the 
same velocity.) For fluids having very high thermal conductivities , 
such as liquid metals, however, heat is probably conducted to or from 
the particle as it moves along its path. This should also occur for 
fluids' having low conductivities if the turbulence is low. In general, 
this effect is to be expected if the ratio of eddy diffusivity to thermal 
diffusivity is small. Because the Peclet number u^D/ (k/pgCp) is a 

measure of the ratio of the average eddy diffusivity at a cross section 
of the tube to the thermal diffusivity, this effect should occur at low 
Peclet numbers. (All symbols sure defined in the appendix.) . In the 
following analysis an attempt is made to account for the heat con du cted 
to or from a turbulent particle as it moves along its path. An analysis 
using somewhat the same general approach is given in reference 7, but 
the assumptions and re stilts differ from those of the present analysis . 

The assumptions used in reference 7 required the eddy diffuslvities for 
heat and momentum transfer to he equal at a Prandtl number of 1, whereas 
experimental evide n ce indicates that this is not always the case. In 
the present analysis , the eddy diffusivity for heat transfer is lower 
than the eddy diffusivity for momentum transfer at low Peclet numbers 
for a Prandtl number of 1. 


Modified Mixing-Length Theory of Heat Transfer 

In the following analysis, eddies are assumed to move transversely 
in the tube between concentric cylindrical surfaces 1 and 2 which are 
separated by the small distance l , the average mixing length at a 
partic ular radius in the tube. Eddies also move longitudinally but these 
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should not affect the radial heat transfer. The eddies originate from an 
instability at either surface 1 or 2 and travel to the other surface -where 
they mix with the fluid. The m-tv-irig might he caused by collisons between 
eddies . 

Equation for turbulent heat transfer . - The case of heat transfer 
from the tube wall to the fluid is considered herein. An eddy of fluid 
originates at the outer surface 1 where it has the time average temperature 
of the fluid at surface 1 t-j_ and arrives at surface 2 with the lower 

temperature tg ' where it mixes with the fluid and assumes the time 
average temperature of the fluid at surface 2 tg • In general, tg ' 
differs from t]_ because of beau conducted out of the eddy between 

surfaces 1 and 2 . The heat transported from surface 1 to surface 2 by 
the eddy is equal to the mass of the eddy multiplied by c g(t„ ’ -t g ) . 

-fir ^ “ 

The turbulent beat transfer per unit area per unit time from surface 1 
to surface 2 is then given by 

<It " PgCpMtg’-tg) (1) 


where f is the fraction of surface 2 on which the fluid in the eddies 
makes contact and v is the average velocity of this fluid toward sur- 
face 2. Equation (l) can be written as 


q t = pgCpfvaC^-tg) 


where 


a s 



( 2 ) 


( 3 ) 
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Inasmuch as l is assumed to be small, the temperature gradient dt/dy 
is approximately constant over l and t-j_ - t 2 can be replaced by 

-1 dt/dy. Equation (2) then becomes 

at “ -PgCp(fvl) a ^ (4) 

or 

Ofc = -P6°p ca (5) 

where c is the coefficient of . eddy diffusivity, the value for which 
depends on the amount and kind of turbulent mixi ng at a point (Ca is 
sometimes replaced by the eddy diffusivity for heat transfer) . 

Equation for turbulent, shear stress . - With the use of the usual 
assumption of-the mixing -length theory, that is , the component of the 
velocity of an eddy in the direction along the axis of the tube does not 
change as the eddy moves between surfaces 1 and 2, the following equation 
for idle turbulent shear stress is obtained: 

T t - p(frl) § (6) 


or 



( 7 ) 


where the eddy diffusivities C in equations (5) and (7) are identical. 
The assumption that they can be considered equal for fully developed flow 
in tubes is further supported by the agreement between the analysis and 
the experiment in the range of Peelet numbers where a equals 1 
(reference 5): 

Determination of a . - In order to obtain a relation for a in 
equation (5) as defined by equation (3), the value for t 2 ' must be . 

determined. Let t ' be the temperature of a point on the eddy at some 
distance s from surface 1. 


The heat transfer per unit area from the eddy at position s is 
equal to h'(t' - t) where h' is the coefficient of heat transfer 
between the eddy and the surrounding fluid. Inasmuch as the dimensions 
and velocity of the eddy are small, the Reynolds number is small and 
laminar heat transfer is assumed. Therefore, . h' is proportional to 
k/&, where 5 is the film thickness or boundary layer associated with 
the heat transfer from the eddy. Then the heat transferred from an eddy 
at s in length ds is kp ds(t* - t)/5 where p is the perimeter 
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of tlie eddy. This heat transfer decreases the temperature of the eddy in 
the direction s and can also he written as -pgCpVA* (dt '/ds)ds, where 

A T is the cross -sectional area of the eddy perpendicular to the 
direction s . Equating these two expressions for heat transfer gives 


3*1 _ ^ / 4 - , _ t i 

flp “ non ffilR U ' 


pgCpVA’S 

As before, dt/dy is considered constant over Z so that 

t = t-t + s 3* 

1 ay 

Substitution of equation (9) in equation (8) and rearrangement gives 
a* 1 . ( feP \ t . ( fcP dt A 

as VpgOpVA'&Z ^pgCpVA’Sy^ 1 dy ) 


(S) 


(9) 


( 10 ) 


Equation (10) can he written as 


3*1 + — £ t « 

ds A'Pe 1 


A'Pe ' (f 1 + dy S ) 


(H) 


where Pe' is the Pec let number for the eddy. Equation (ll) is a first 
order linear differential equation In the variables t * and s The 
solution of equation (ll) with the initial condition that t ' = t]_ 
when s * 0 is 


dt A*Pe* dt 

ay 


t . 10.1/ a x w u.u . _ 

* = S -5T1 — + + — sr: e 


A'Pe* dt 

ay 


"G-'Pe’) 5 


ay p ay ' ' p 

Eor s .as 2, t* = t 2 *, and also t^_ - t 2 = - 2 dt/dy, then 





a = 


A'Pe 1 

2p 



( 12 ) 


(13) 


In order to use equation (13) , it is necessary to mate some 
assumptions concerning the quantities p. A’, and Pe*. The first 
assumption made is that the dimensions of the eddy are proportional to 
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the mixing length l, or A' is proportional to p 2 which is pro- 
portional to l 2 . An assumption must also he made for Pe*. As a 
first approximation, Pe' is assumed proportional to Pe, the Peclet num- 
ber for the tube, and the results will be Compared with experimental 
results. Equation (13) becomes, by introduction of these assumptions. 


a = b Pe 



(14) 


where b is a constant to be determined experimentally. Equation (14) 
gives the desired relation for a to be substituted in equation (5) . 


Calculation of Nusselt Number , Reynolds Numbers, and Peclet Number 

for Uniform Heat Transfer and Uniform Fluid Properties 

When the molecular shear stress is added to the turbulent shear 
stress given by equation (7) and the molecular heat transfer is added 
to the turbulent heat transfer given by equation (5), the following 
equations for shear stress and heat transfer per unit area are obtained: 

t = U. + pe) 0 (15) 

1 - - (k + pgOp<tt) jjp (16) 

where a is given by equation (14) . The eddy diffusivity C is obtained 
from the following expressions which are experimentally verified in refer- 
ences 5 and 10: 


C = n 2 uy (17) 

for flow close to the wall (y+<26) and the Karmen relation 

£=K 2j SiZSdfj (18) 

(d 2 u/dy 2 ) 

for flow at a distance from the wall (y + >26), where n and tc are 
experimental constants having the values 0.109 end. 0.36, respectively. 

By integration of equation (15) for the region at a distance from the 
wall with the viscous shear stress neglected, von Kahnan obtained, for 
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•uniform properties, the following linens ionless equation (reference 10): 

(19) 


U + = i 




-H C 


Equation (18) becomes , on substitution of the first and second derivatives 
from equation (19) , 




( 20 ) 


for flow at a distance from the wall. For flow close to the wall, equa- 
tion (17) can be written in dimensionless form as 



n2 u+y+ 


( 21 ) 


where the values of u + are obtained from, the generalized velocity dis- 
tributions (constant fluid properties) given in reference. 10. Values 
of e from equations (20) and (21) are substituted into equation (16) 
for obtaining temperature distributions. Equation (16) can be written 
in dimensionless form as 


_SL = . 

90 


1 



( 22 ) 


The dimensionless quantities are defined in the list of symbols. Equa- 
tion (22) can be written in integral form for calculating temperature 
distributions as 


n y 


v - t " - 


iAsq 


i + ^ ac' 


dy* 


(23) 


where 


C = 2 g 
r 0 (n/p) Re 


(24) 


The expression for Re will be given in the section Reynolds and 

Peclet numbers. An equation for q/qo will be obtained in the following 

section. 
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Variation of q/q^ across tube. - Die variation of heat transfer 

per unit area across the tube can be obtained by taking a heat balance 
on an annulus of fluid having differential length dx, outside radius r, 
and thickness dr. The heat added to the annulus must be equal to the 
heat transferred into the annulus minus the heat transferred out of the 
annulus, or 

r? 4* 

2itr dr pguc^ ^ = 2jtrq - 2ic(r - dr) (q - dq) 


or, neglecting higher order differentials, 

d(rq) = pgucp || rdr (25) 

If dt/dx is constant across the tube, as must be the case for a fully 
developed temperature distribution with uniform heat transfer, the 
following heat balance can be written for the cross section as a whole: 

*r 0 2 Pg%c p H = 2 I tr 0 q 0 (26) 

Dividing equation (25) by equation (26) gives 


Ub r 0 2(30 


(27) 


Replacing r by (r Q - y) and writing equation (27) in dimensionless 
form yields 


d 


- f + > £ 


+ (r 0 + - y + ) + 

= - 4u + — — dy + 


Re 


(28) 


Integrating gives 

p(ro + - y + )q/<io 


vJr, 


rr 


0 


- y+) - 


4 

Re 


u + (r Q + - y + ) dy + 


Oo 


or 


r 0 


90 r 0 + - y + Re(r 0 + - y + ) 


r y 


Oo 


u + (r 0 + - y + ) ay+ 


(29) 
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Equation {29) gives the desired relation 'Which is to he substituted into 
equation (23) for the variation of q/q q across the tube. 

Nusselt numb ere . - The bulk temperature t^ for constant fluid 
properties is defined by 


% 


.1 


tur dr 


J : 


(30) 


ur dr 


Equation (30) can be written in dimensionless form as 

’) u + (l - y , ) a y I 


j- n j. « 

t o ' \ " 


Jo (t °” " tn) 


Jo* 


(31) 


+ (i - y')<3y’ 


where tQ n - t" is obtained from equation (23 ) , but 

*0 


h s 


*0 - V 


therefore 


U US M = 2 

T V - V 


(32) 


(33) 


Reynolds and Peclet numbers . - The Reynolds number Re is given by 

(34) 


Re - = 2u b + r 0 + 


where u^ is given in reference 5 as 


rtf*. 


V = 


(r n + ) 


Wo 


+ f + +N, + 

« (r 0 - y )<3y 


(35) 
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The Peclet number is calculated from 

Pe = (Re) (Pr) (36) 

Calculation of Husselt number as a function of Reynolds or Peclet 
number and Prandtl number . - In the calculation of Nusselt number, the 
values for r Q + ‘ and Pr are to he assumed first . With the use of the 
u + against jr - curve from reference 10 (fig. l) , can he calculated 

from equation (35) . Then Re can he calculated from equation (34) and 
Pe from equation (36). Hie next step is to obtain e/(p/p) as a function 

of y + from equations (20) and (21) , and €’ from equation (24). Then, 
q/q^ is obtained from equation (29) and y' from y' = y+/rQ+ for each 

value of y+. Hie quantity a can be calculated from equation (14) . The 
generalized temperature difference to" - t" can then he calculated as 
a function of y' from equation (23). With the temperature distribution 
known, to" - t^" can he obtained from equation (31) and the Husselt num- 
ber from equation (33) . 


RESULTS AND DISCUSSION 

Variation of heat transfer per unit area across tube for uniform heat 
transfer at surface . - The variation of heat transfer per unit area with 

y + for various values of Vq + is shown in figure 2 . The curves indicate 

that the heat transfer per unit area varies nonlinear ly from a maximum 
at a point near the wall to zero at the tube center. For low values of 

Tq + , the maximum heat-transfer per uni t area is attained at a point 
slightly removed from the wall rather than at the wall. This is because 
the heat -transfer area decreases as the distance from the wall increases 
and thus, even though the total heat transfer decreases as the distance 
from the wall increases, the heat transfer per unit area increases 
slightly for some conditions . 

Experimental determination of constant b; Nusselt numbers for gases . - 
In order to calculate relations among Nusselt number, Reynolds or 
Peclet number, and Prandtl number it is necessary to determine exper- 
imentally the constant b in equation (14=) . This constant is obtained 
from experimental data for heat transfer to air. in reference 5 and plotted 
in figure 3. Inasmuch as the data were obtained at high ratios of wall 
to fluid bulk temperature, the temperature at which the fluid properties 
in the Nusselt and Reynolds numbers are evaluated is important, mhg data 
indicate that, for Reynolds numbers below lt>,000, the best correlation 
was obtained by evaluating the properties at the bulk temperature j for 
Reynolds numbers between 15,000 and 30,000, equally good correlation was 
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was obtained by evaluating the properties at the bulk temperature or at 
t 0 . 4 "Where 1 0>4 _ = 0.4(t 0 t^) + t^,; and for Reynolds numbers above 

30,000, the best correlation was obtained by evaluating the properties 
at tQ^. In. figure o the properties are evaluated at the hulk tem- 
perature for Reynolds numbers below 30,000 and at t^ ^ for Reynolds 

numbers above 30,000. A justification for evaluating the properties at 
the b ulk temperature rather than at "to. 4 for the lower Reynolds numbers 

is found in reference 11, in which it is shown that for la minar flow the 
properties should be evaluated at a temperature between the bulk tem- 
perature and the_ temperature at the center of the tube . It would then 
be expected that in the -transition region the properties should be 
evaluated at a temperature between this temperature and tQ 

Predicted curves from the analysis are included in figure 3. It is 
apparent that for a = 1 . (dashed curve) , the predicted Husselt numbers 
at low Reynolds or Peclet numbers are too high, whereas the solid curve 
for which a is calculated from equation (14) is in good agreement with 
the experimental results. The value of the constant b in equation (14) 
is 0.000153 as determined from the experimental data. 

In order for the analysis to apply, the Reynolds number should be 
well above the Reynolds numbers at which transition from laminar to tur- 
bulent flow takes place. The velocity-distribution data for air (fig. l) 
indicate no effebt of transition for Reynolds numbers as low as 8000. On 
the other hand, -the temperature -distribution data for air from refer- 
ence 5 indicate a considerable effect at that Reynolds number so that the 
data in figure 3 at Reynolds numbers in the vicinity of 10,000 should be 
suitable for determining the constant a. 

Liquid-metal heat transfer . - Predicted curves for Husselt number 
against Peclet number for a Prandtl number of 0.01 are shown in figure 4. 
The curve for a = 1 is substantially the same as those obtained in 
references 1 and 2. The solid curve was obtained by computing a from 
equation (14) where the value of b (0.000153) was obtained from the data 
for flow of air in figure 2. Curves representing data for mercury (data 
obtained from. Musser and Page, reference 2) and lead bismuth from refer- 
ence 3, are also shown. It Is shown in reference 2 that for Prandtl num- 
bers below O’.! the curve of Nusselt number against Peclet number Is 
practic al Ly independent of Prandtl number so that the predicted curves 
for Pr = 0.01 should be comparable to the data for mercury and lead ' 
bismuth. It is apparent from figure 4 there is improvement in the 
agreement between experimentally and analytically determined values ' 
when a Is computed from equation (14) although there Is considerable 
scatter in the data. Data at higher Peclet numbers, in the region where 
the two predicted curves begin to merge, would be of value in order to 
establish further the validity of the analysis . This data would indicate 
whether or not a actually approaches 1 at high Peclet number for liquid 
metals. 
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An empirical equation of the type given in reference 2 is also 
plotted in figure 4. This equation is 


Nu = 6.3 + 0.000222 (Pe) 1 * 3 (37) 

Equation (37) represents the predicted line with a computed from equa- 
tion (14:) with good accuracy for the range of Peelet numbers shown in 
figure 4. . 


SUMMARY OF RESULTS 

The following results were obtained from the analytical investigation 
of turbulent heat transfer at low Peelet numbers: 

1. The mixing -length theory was modified in order to include the 
heat conducted to or from an eddy as it moves transversely. In order to 
utilize the analysis, a constant was deter mine d experimentally. 

2. When the experimental constant was determined from a point on the 
curve of Nusselt number for gases, the analytical results were brought into 
better agreement with the data for both gases and liquid metals at low 
Peelet numbers. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX - SYMBOLS 


cross -sectional area of eddy perpendicular to radial direction, ft^ 
experimental constant > 

constant of integration 

specific teat at constant pressure, Btu/lb (°P) 
inside diameter of tube, ft 

fraction of surface 1 or 2 on which eddies are arriving from the 
other surface 

acceleration due to gravity, 52.2 ft/sec 2 

coefficient of heat transfer betveen wall and fluid, q^/ (t Q -t b ) , 
Btu/(sec)(ft 2 )(°P) 


coefficient of heat transfer between eddy and fluid, Btu/(sec) (ft^) (°p) 

thermal conductivity of fluid, (Btu) (ft)/(sec) (ft 2 ) (°P) 

average mixing length at a particular radius, ft 

Nusselt number, ho/h 

experimental constant 

Peclet number for tube, (pgCp)j 

Peclet number for eddy, v5/^k/ (pgCp)j 
Prandtl number, Cpgp./k 
perimeter of eddy, ft 

total rate of heat transfer per unit area toward tube center, 

(Btu/ (sec) (ft^) 


rate of heat tr ans fer per unit area toward tube center by turbulent 
motion, Btu/(sec) (ft 2 ) 

total rate of heat transfer at wall toward tube center per unit 
area, Btu/(sec) (ft^) 

Reynolds number, pu^D/u 
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r radius, distance from tube axis, ft 
Tq inside tube radius, ft 

Tq 4- tube radius parameter. 


m/p r ° 


S 


t 


t 


0 


0.4 


t 


2 


t* 


V 

% 

t" 

V 

V 




V 

v 

x 

y 


distance toward tube center from surface 1 where eddy originates, ft 
time average temperature at a point, °F 
wall temperature, °F 

0.4 (t 0 -t„) + t,, 

time average temperature .at surface 1, where eddy originates, °F 

time average temperature at surface 2, where eddy mixes with 
fluid, °F 

temperature of a point on eddy a distance s from surface 1, °F 
temperature of eddy at surface 2, °3? 

bulk or average temperature at a cross section of tube, °P 
dimensionless temperature parameter, (k/q q rg)t 
dimensionless wall temperature parameter, (k/q^ rQ)tQ 
dimensionless bulk tenperature parameter, (k/q^ Tq)^ 
time average axial velocity at a point in fluid, ft/sec 
bulk or average velocity at a cross section in tube, ft/sec 
velocity parameter, u /"\/ T o/f> 

bulk velocity parameter, '%/''/ Tq/p 
radial velocity of eddy, ft/sec 
distance from tube entrance, ft 
distance from tube wall, ft 

wall distance parameter. 



y 
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y* y/r 0 

a (V “ - V' € / c li 

5 film thickness or thickness of boundary layer on eddy, ft 

€ coefficient of eddy difrusivity for momentum, ft 2 /sec 

coefficient of eddy diffusivity for heat transfer, ft 2 /sec 

«/(«b r 0 ) 

K Karmen constant 

H absolute viscosity of fluid, (lb) (sec)/ft 2 

p mass density, (lb) (sec 2 )/fu 4 

T shear stress in fluid, lb /ft 2 

T^_ turbulent shear stress in fluid, lb /ft 2 

Tq shear stress in fluid at ■wall, lb /ft 2 
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Figure 5. - Variation of Suseelt number with Reynolds number for heat transfer to air. Prandtl number, 0.73. Physical pro- iP 
perties for data evaluated at bulk ten£>erature for Reynolds number below 30,000 and at temperature, t n a * 0,4 (t^-tv) -f tuJ* 0 
for Reynolds number above 30,000, * 1 
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